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The thermally induced acoustic effect occurring in a regenerator with a temperature 
gradient has been observed on a particularly built regenerator experimental set-up 
attached to a Stirling cryocooler system. To describe this thermoacoustic phenom¬ 
enon, a parity simulation network model has been established in the form of an oscil¬ 
lation circuit with positive feedback amplification. The model vividly shows the influ¬ 
ence on the whole system of the change in flow impedance in the regenerator caused 
by the change in temperature gradient. 
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Nomenclature 



A 

Cross-sectional area (m 2 ) 

Superscripts 

C 

Volume capacitance of fluid (m 5 N~') 

> 

Composite 

C d 

Motion damping factor of displacer (N s m ') 


Derivative of first order 

F 

Amplification ratio of feedback circuit 


Derivative of second order 

G 

Voltage ratio of secondary in transformer 

—► 

Phasor 

/ 

Equivalent current in network 



K 

Rigidity of spring of displacer (Nm H ) 

Subscripts 

L 

Equivalent inductance of fluid (N s 2 m -5 ) 

M 

Mass of displacer (kg) 

ab, be 

Two parts of secondary in voltage 

P 

Pressure wave (N m 2 ) 


transformer 

Q 

Flow flux (m 3 s~‘) 

e, w, d 

Cold space, warm space and displacer of 

R 

Equivalent resistance of fluid (N s m~ 5 ) 


cryocooler 


Gas constant (m 2 s 2 K -1 ) 

f 

Feedback 

S 

Complex frequency for Laplace 

g 

Pneumatic space in cryocooler 


transformation 

i, o 

Inlet and outlet 

T 

Temperature (K) 

m 

Matrix of regenerator 

U 

Equivalent voltage in network 

Volume (m 3 ) 

r 

Regenerator 

V 

rod 

Driving rod of displacer 

X 

Displacement of displacer (m) 



z 

Equivalent impedance of regenerator 



Greek letters 



P 

Amplification gain 



A T r 

Longitudinal temperature difference in 
regenerator (K) 



e 

Filling factor of matrix 



OJ 

Angular frequency 
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From a thermoacoustic point of view, it is now being reco¬ 
gnized more and more commonly that in a regenerative heat 
engine the conversion between heat and work is essentially 
carried out by the coupling of a longitudinal pressure wave 
propagating in the working fluid with a longitudinal tem¬ 
perature gradient of the matrix via transverse heat transfer 
in the regenerator, i.e. thermoacoustic conversion by the 
regenerator. This causes, in the case of a positive tempera¬ 
ture gradient along the line of pressure propagation, the 
amplification of acoustic energy, while it generates heat 
pumping in the case of a negative temperature gradient '~ 3 . 
In more detail, the thermoacoustic effect in a regenerator 
is described as the increment of acoustic energy in the 
direction of the temperature gradient due to heat addition 
by a transverse temperature wave, modulated by gas parcel 
displacement along with the propagation of a pressure 
wave 4 . After successful modelling of pressure wave propa¬ 
gation in an isothermal regenerator by a passive network 
of fluid flow reactance (see Figure la) 5 , recently the ther¬ 
moacoustic conversion in a regenerator has been charac¬ 
terized by an active network made up of acoustic imped¬ 
ance and a source element reflecting the contribution of gas 
parcel velocity variation during its reciprocating displace¬ 
ment along the longitudinal temperature gradient (see Fig¬ 
ure lb) 6 . 

It has also been found experimentally that, in a split cycle 
Stirling cryocooler, the cold space temperature significantly 
influences the working process of the cryocooler system 7,8 . 
From the thcrmoacoustic point of view, this effect reflects 
the interaction between the longitudinal temperature gradi¬ 
ent of the matrix and wave propagation in the working gas. 
A detailed investigation would reveal more about the mech¬ 
anism of conversion between heat and work in the regener¬ 
ator. For this purpose, a series of experiments have been 
carried out on a specially designed split cycle Stirling cryo¬ 
cooler to investigate the mechanism of the conversion. 
Based on network modelling of the regenerator, the idea of 
parity simulation proposed first at MIT 9 has been developed 
to apply to a parity simulation circuit for a regenerator. The 
endeavour to combine the thermoacoustic viewpoint with 
the fluid network modelling theory 10 to provide a more con- 



b 

Figure 1 Equivalent network of wave propagation for regener¬ 
ator. (a) Passive form for isothermal system; (b) active form for 
regenerator with temperature gradient 



Figure 2 Schematic diagram of experimental set-up 

venient engineering tool is also one of the objectives of 
this paper. 

Experimental set-up and main results 

An experimental split Stirling cryocooler (shown in Figure 
2) with a controllable displacer (driven by a linear motor) 
and an adjustable longitudinal temperature difference ATj. 
exerted between the ends of the regenerator, has been built 
specially for this investigation. The temperature in cold 
space Tj, can be determined by a liquid nitrogen source and 
an electric heater both pressed tightly against the head cup 
of the cold space. The linear motor attached to the displacer 
is fed with square wave signals by a frequency regulator. 
During the experiments, by controlling the temperature T c 
and driving the displacer to move in a sinusoidal pattern, 
the pressure waves in the cold and warm spaces P c , P w and 
the displacement of the displacer X d , have been measured 
accurately under various temperature differences A T r . All 
these steps have been carried out under the condition that 
the compressor is out of operation. Two examples of the 
measured pressure and displacement waves with a selected 
temperature difference A T r , after the noise being filtered, 
are shown in Figures 3 and 4, while all the results are sum¬ 
marized in Figures 5 and 6. The experimental results reveal 
the following main facts. 


Xj(mm) P(Kgf/cm 2 ) 



Figure 3 Results of X d , P E and P w fluctuation (after filtration) 
measured at AT R = 0 K (7 E = T w = 300 K) 
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Figure4 Results of X d , Pe and P w fluctuation (after filtration) 
measured at A T p = 163 K (T E = 122 K, T w = 285 K) 


Degrees Degrees 



At„«o 


Figure 5 Phase shift between X d , P E and Pvv terms of A T„. 
o, P E A P w ; o, P w A X d ; x, P E A X d 


Ap E (Kgf/cm 2 ) A Pw (Kgf/cm 2 ) 



AW) 


Figure 6 Amplitude variation of P E and PM with AT b . o, A P j o, 
Ap m (charging pressure = 11.0 kgf cm -2 ) 

The longitudinal temperature difference A T r imposes a 
significant influence on the pressure wave propagation 
induced by the motion of the displacer in the regenerator, 
which can be characterized not only by phase shifts 
between P e , P w and X d , but also by the increment of the 
amplitudes of P e and P w . Specifically, at zero A T r , P e has 
a phase lead over P w of =180° (shown in Figure 3), and 
the amplitude of P w is very small compared to P e because 


the warm space is connected with the compression space 
by a linking pipe and appears to be almost the full portion 
of the dead volume of the cryocooler system, while the 
volume of the cold space occupies only a very small portion 
of it. With the appearance of AT r , however, the phase lead 
of P e over P w drops rapidly, and P e goes ahead of P w by 
=90° even if the A T r is small (shown in Figure 5). During 
the process of A T, changing from 0 to 110K, the phase 
difference between P e and P w reduces slowly from =90° 
to 60°, whereas the amplitudes of P w and P e (especially 
P w ) obviously increase (shown in Figures 5 and 6). When 
the temperature difference A T t attains the level of 120 K, 
the phase of P w shifts abruptly and comes into phase with 
P e . To compensate the violent phase shifting of P w , its 
amplitude reduces. 

As the temperature difference A7' r continues to increase, 
the in-phase relationship between P e and P w does not obvi¬ 
ously change and the amplitudes of these two pressure 
waves rise together. Thus there are two points where abrupt 
phase shifting between P e , P w and X d (shown in Figure 5) 
occurs with a corresponding variation in the amplitudes of 
P e and P w (shown in Figure 6) in the experimental tempera¬ 
ture difference A T, range; and here the second point is more 
noteworthy. In contrast to the changing phase relation 
between P e and P w , the phase relation between P e and X d 
shows no obvious variation and P c continues to have a 
phase lead over X d of =90° in the whole experimental 
range. 

Parity simulation based on network 
modelling of regenerator 

To simulate the wave process caused by the motion of the 
regenerator, the displacer’s cylinder and the movable regen¬ 
erator inside it can be described as a cylinder with a recipro¬ 
cating piston in it, which divides the cylinder into two 
unequal separated spaces, V w and V e ( V w A V e ), both of 
them being connected by a bypass tube with a flow imped¬ 
ance equivalent to that of the regenerator. A closed-core 
transformer with a specific point tapped secondary can be 
used to model the pressure wave effect induced by the 
motion of the piston, and the displacement rate of the 
regenerator X d is modelled as the current of the trans¬ 
former’s primary / xd which provides two voltages (7 ab and 
U cb in the secondary, each with the opposite phase to the 
other to simulate the pressures P w and P c in the case of 
zero A T v . According to the already known network 
expression for wave propagatibn in a stationary regener¬ 
ator 6 , a parity circuit of the pressure wave effect produced 
by the motion of the regenerator can be illustrated as in 
Figure 7a, in which the tap point of the secondary of the 
transformer shifts to the end near C w to give an amplitude 
of U ab much larger than that of U cb . Here, the capacitors 
C c , C r and C w represent the mass accumulation capacity of 
the cold space, regenerator and warm space, respectively, 
and the inductor L, and resistor R r correspond to the effect 
of inertia and viscosity of the working fluid in the regener¬ 
ator. The circuit gives the following relations 


U* = 

U cc + u r 

(1) 

U cb = 

U c „ + U T 

(2) 

^ce 

/ r + / cw 

(3) 


which agree with the conservation principle in a regenerator 
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Figure 7 Parity simulation network for isothermal regenerator, 
(a) Circuit of network; (b) graph for vector analysis 


and two adjacent spaces. In the above equations, voltage 
U, represents the pressure drop across the regenerator, and 
the currents 7 CC , / r and / cw stand for the mass flow rate of 
gas in the three elements. Since U ah is in phase with / xd , in 
contrast U cb , U cc and l/ cw are nearly out of phase depending 
on U r in Figure 7b. Quantitatively, by applying the prin¬ 
ciple of equivalent transformation of a two port active net¬ 
work to the calculation, the transfer function between U ce 
and U cw can be obtained as 


U ce (s) _ G + C w x Z, x s 
U c Js ) = (1 -G) + C e xZ r xs 

here Z r is the total impedance of the regenerator 

R. + L r X s 

Z _- i -1- 

r (l - or x C T x L r ) + R r x L r x s 

and G is determined as 


(4) 


(5) 


_ |£4bl 

It/J 


( 6 ) 


from which it can be seen that the phase relation between 
P e and P„ is significantly determined by the impedance Z r 
of the regenerator. In accordance with the fact that under 
A T r = 0 K the flow impedance of the regenerator is domi¬ 
nated by the fluid capacitance Cj. (reference 5), so U CT simu¬ 
lating the pressure drop component contributed by this 
capacitance should be larger than t/ Lr caused by fluid iner¬ 
tia. Phasor analysis of the phase relation between U„., U CVI 
and / xd in Figure 7b shows good agreement with the exper¬ 
imental situation where A7’, = 0 K (seen in Figure 3). 

The appearance of a temperature gradient on the regener¬ 
ator makes the network active because a transverse heat 



Figure 8 Parity simulation network for regenerator with tern 
perature gradient, (a) Circuit of network with oscillation; (b) 
graph for vector analysis 


flux is modulated by the pressure wave and this contributes 
to sound generation, i.e. pressure amplification under the 
condition of a positive temperature gradient. This coincides 
with the situation in the experiments, so it is clearly shown 
that, when the temperature gradient is raised, the ampli¬ 
tudes of P e and P w increase separately in varying degrees, 
while the disturbance X d almost does not change as shown 
in Figures 3 and 4. This phenomenon is analogous to the 
situation in a oscillation circuit with a positive feedback 
amplifier, which modulates a d.c. voltage to an a.c. one, 
with its input being amplified at the same time. To describe 
this effect, the source in the network shown in Figure lb 
is simulated by the positive feedback amplifier circuit 
shown in Figure 8a, in which the d.c. power source F s 
corresponds to the heat input and its power level increases 
along with augmentation of the temperature gradient. The 
voltage drop U cc on C e is amplified by a two-stage amplifier 
whose output is fed to the parallel circuits B-O and C-O. 
A couple of transistors are needed to satisfy the phase 
relation for the positive feedback amplifier circuit because 
the output of a single transistor is in opposite phase to its 
input. The positive feedback loop is composed of the paral¬ 
lel circuit and part of the secondary of transformer A-B 
(shown in Figure 9), so that the output of the amplifier 
circuit is fed back to the input. From analysis of the positive 
feedback loop 


L, Z' 


n 




_ 

Ur 

= C E 

1 

L-2 \ 

L Zr - 

U 0 

(P E ) 


C w = 

= j 

(PJ 







Figures Positive feedback loop of circuit for regenerator 
with oscillation 
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F 


U f L, 
~=1 +^xs 


(7) 


can be obtained, in which impedance Z' includes the com¬ 
bination of Z r , C w and L 2 , while L, and L 2 are the inductance 
of the two parts of the secondary A-B, B-C. In an ideal 
case, L,-L 2 = 0 , so F= 1. If the temperature difference 
across the regenerator is initially small, the gain of the 
amplifier circuit is strong enough to contend with the effect 
of the secondary of the transformer, and so the voltage 
relations (1) and (2) must be changed to (8) and (9) 

U ce = 11.^+U r (8) 

£4 = t7 bc + (7 CW (9) 

which means that the amplifier circuit becomes a determin¬ 
ing component in the power supply of the system. The 
phasor analysis shown in Figure 8b gives the phase relation 
between £/ ce , £/ cw and / xd , which is in good agreement with 
the case of 0 < at, < 120 K. In the process where A T t 
increases up to about 120 K, I/ cc continues to be amplified 
to keep C/ cw amplified, while the phase shift between U ce 
and U cv/ is slowly decreased. The reason behind this change 
will be given in a later analysis. Evidently, the increase of 
A T, corresponds to the increment of heat input, i.e. aug¬ 
mentation of the mains voltage V s in the feedback circuit. 

Next, we go further into the situation where A T r becomes 
larger than 120 K. In the network in Figure 8a, the flow 
impedance of the regenerator Z r would change because the 
flow resistance of the matrix decreases continuously along 
with the cold space temperature reduction, while the flow 
inductance is steadily augmented and the fluid capacitance 
remains almost unchanged (R, * /a, L, oc p , C r * 
l/(/? g x P) = constant). When AT r reaches —120 K, the 
flow inductance of the regenerator increases so as to be 
large enough to compensate thoroughly the effect of fluid 
capacitance; a kind of parallel resonance occurs in the 
regenerator which shows a relatively pure resistance 
characteristic. Thus the phase relation between (7 ce and U CK 
undergoes a sudden change, so that U cw is in phase with 
U cc according to the phasor analysis shown in Figure 10. 
So it can be concluded that, along with the establishment 
of a temperature gradient, the flow impedance of the regen¬ 
erator varies gradually from the capacitance-dominated 
state to the resistance-dominated one. In the process where 
AT, continues to increase, the resonance condition remains 
and the amplification effect steadily grows while the heat 
input increases. It is regrettable that the authors failed to 
complete the experiment reported above to attain a cold 
space temperature of 77 K. But from the results in another 


O C Ubc b U ba a 


Uoc <PJ ( 

UcB 

UpB^r) .. 

U CA 

(U CW ) 







(U ce ) (P E )U OA 


Figure 10 Vector analysis graph for regenerator network at its 
resonant state 
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Figure 11 Results of P w and X d fluctuation measured at 
T e = 94 K 


similar experiment 8 shown in Figure 11, it may be deduced 
that the departure from the resonance state would occur 
in the case of T e = 94 K, corresponding to AT r = 200 K, by 
comparing previous values with the results in the range 
0 < AT r < 120 K. 

Thus, a parity simulation circuit based on the thermo¬ 
acoustic network modelling of a Stirling cryocooler regen¬ 
erator has been established and shown to be in good agree¬ 
ment with the experimentally revealed effect of temperature 
gradient on pressure wave propagation. 

Combining the analysis above with the results of another 
experiment reported in reference 8, in which the variation 
of the damping factor C d of the displacer with respect to 
refrigeration temperature was determined by an on-line 
method 7 , i.e. measurement in a working Stirling cryocooler 
system, it can be seen that the dependence of the flow 
impedance of a regenerator on the cold space temperature 
is significantly relevant to the working performance of the 
cryocooler. One especially important result of the exper¬ 
iment in reference 8, the change in shape of the damping 
factor C,], expressed in the motion equation (10), for the 
displacer as the refrigeration temperature was dropping 

MxX;'+C d xX d + A-xX d = (P g - P w ) x A rod (10) 

+ (P w -T c )x(A m xe + A d -A m ) 

is shown in Figure 12. Obviously, the change in motion 
damping with cooling down has strong non-linear behav- 



80 100 120 140 160 180 200 220 240 260 280 300 


T E (K) 

Figure 12 Changing pattern of C d with T E 
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iour, shown as the stepping down of damping when the 
refrigeration temperature is higher than 190 K, then the pro¬ 
gressively accelerating decline until the minimum is 
reached at =130 K and, finally, the sudden step rise when 
the refrigeration temperature is below 120 K. 

Based on the experimental facts presented in this paper 
and the above discussion, it is easy to understand the reason 
behind the singular change of damping. In fact, the motion 
damping changes along with the dropping refrigeration 
temperature in two opposite ways: one trend is to decrease, 
due to the monotonic decrease of gas viscosity, being rep¬ 
resented as the stepwise reduction of flow resistance of the 
matrix in the regenerator; while another trend is to increase 
on the appearance of the activeness of the regenerator, 
especially when the resonance state of the regenerator 
occurs. Evidently, at a refrigeration temperature lower than 
160K(Ar r >120 K), the motion damping of the displacer 
is significantly determined by the state of flow impedance 
of the regenerator. From that, it can be seen that the 
activeness induced by the thermoacoustic effect in a regen¬ 
erator is a key link between the so-called negative feedback 
relationship of the cold space temperature and the perform¬ 
ance of the cryocooler, proposed first in reference 4. 


Conclusions 

From the experimental results and the above analysis, a 
thermally induced acoustic effect in a regenerator with a 
temperature gradient, namely the activeness of the regener¬ 
ator, has been found and successfully simulated by a parity 
simulation circuit with positive feedback amplification. The 
activeness of the regenerator can be generated not only by 
the modulation of pressure wave propagation, e.g. in a 
stationary regenerator, but also by the motion of the regen¬ 
erator, e.g. the regenerator in a Stirling cryocooler. More¬ 
over, it would change along with the temperature gradient 


in the form of variation of flow impedance of the regener¬ 
ator, as observed in our experimental range where the flow 
impedance of the regenerator changes from the flow capaci¬ 
tance-dominated state to an important flow resistance-domi¬ 
nated resonance state as the temperature gradient across the 
regenerator steadily increases. As a result, the performance 
of the cryocooler has been greatly affected by this 
activeness. 
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